Observation of two charged bottomonium-like resonances in T(5S') decays 
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We report the observation of two narrow structures in the mass spectra of the 7r*T(n5') (n = 1,2,3) 
and 7r^/ii,(mP) (m = 1, 2) pairs that are produced in association with a single charged pion in T(5S') 
decays. The measured masses and widths of the two structures averaged over the five final states 
are Mi = (10607.2 ± 2.0)MeV/c^ Ti = (18.4 ± 2.4) MeV and M2 = (10652.2 ± 1.5)MeV/c^ 
r2 = (11.5 ± 2.2) MeV. The results are obtained with a 121.4 fb~^ data sample collected with the 
Belle detector in the vicinity of the T(5S') resonance at the KEKB asymmetric-energy e+e~ collider. 

PACS numbers: 14.40. Pq, 13.25. Gv, 12.39.Pn 



Recent studies of heavy quarkonium have produced 
a number of surprises and puzzles [ij, including some 
associated with T(bS) decays to non-BB final states. 
The Belle Collaboration reported the observation of 
anomalously high rates for T(5S') — s- T(??S')7r+7r^ {n — 
1,2,3) and T(55') ^ hh{mP)TT+ ii~ (m = 1,2) Q 
transitions. If the T{nS) signals are attributed entirely 
to T(5S') decays, the measured partial decay widths 
r[T(55) ^ T(n5)7r+7r-] - 0.5 MeV are about two or- 
ders of magnitude larger than typical widths for dipion 
transitions among the four lower T{nS) states. Further- 
more, the processes T{5S) — > /ib(mF)7r+7r~, which re- 
quire a heavy-quark spin flip, are found to have rates 
that are comparable to those for the heavy-quark spin 
conserving transitions T{5S) — > T(nS')7r+7r~ [3|. These 
observations differ from apriori theoretical expectations 
and strongly suggest that exotic mechanisms are con- 
tributing to T{5S) decays. We report results of res- 
onant substructure studies of T(55') — )• T(nS')7r+7r~ 
(n = 1,2,3) and T(55) hb{mP)Tr+7r- (m = 1,2) de- 
cays We use a 121.4 fb""'^ data sample collected on or 
near the peak of the T{5S) resonance (-^/s ^ 10.865 GeV) 
with the Belle detector at the KEKB asymmetric energy 
e+e~ collider 

The Belle detector is a large-solid-angle magnetic spec- 



trometer that consists of a silicon vertex detector, a 
central drift chamber, an array of aerogel threshold 
Cherenkov counters, a barrel-like arrangement of time- 
of-flight scintillation counters, and an electromagnetic 
calorimeter comprised of CsI(Tl) crystals located inside a 
superconducting solenoid that provides a 1.5 T magnetic 
field. An iron flux-return located outside the coil is in- 
strumented to detect mesons and to identify muons. 
The detector is described in detail elsewhere @. 

To reconstruct T(5S') ^ T(7iS')7r+7r- , T{nS) ^1+ 
candidates we select events with four charged tracks 
with zero net charge that are consistent with coming 
from the interaction point. Charged pion and muon 
candidates are required to be positively identified. Ex- 
clusively reconstructed events are selected by the re- 
quirement |Mmiss(7r+7r") - M{^l+^-)\ < 0.2 GeV/c^, 
where Mmiss ("'''' tt") is the missing mass recoiling 
against the tt+tt" system calculated as Mi„iss(7r+7r^) = 

(Ec.m. ~ -p*+^-, Pern, is the center-of-mass 

(cm.) energy and and are the energy 

and momentum of the tt+tt^ system measured in the 
cm. frame. Candidate T(5S') T{nS)Ti^TT~ events 
are selected by requiring |Mi„iss(7r+7r~) — tox(„5)| < 
0.05GeV/c^, where mx(nS) is the mass of an T{nS) 
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state 0- Sideband regions are defined as 0.05GeV/c^ < 
\M^iss{n+TT-) - mrins)\ < O.lOGcV/c^. To re- 
move background due to photon conversions in the 
innermost parts of the Behe detector we require 
A'/2(7r+7r-) > 0.20/0.14/0.10 GeV/c2 for a final state 
with an T(1S'),T(25),T(35'), respectively. 

Amplitude analyses of the three-body T{5S) — > 
T(nS')7r+7r~ decays reported here are performed by 
means of unbinned maximum likelihood fits to two- 
dimensional M^[T{nS)Tr+] vs. AP[T{nS)Tr-] Dalitz dis- 
tributions. The fractions of signal events in the sig- 
nal region are determined from fits to the correspond- 
ing Af„iiss(7r+7r~) spectrum and are found to be 0.937 ± 
0.015(stat.), 0.940 ±0.0Q7(stat.), 0.918 ± 0.010(stat.) for 
final states with T(15'), T(2S'), T(35'), respectively. The 
variation of reconstruction efficiency across the Dalitz 
plot is determined from a GEANT-based MC simula- 
tion Q. The distribution of background events is de- 
termined using events from the T(nS) sidebands and 
found to be uniform (after efficiency correction) across 
the Dalitz plot. 

Dalitz distributions of events in the T(2S') sidebands 
and signal regions are shown in Figs, [ija) and[ljb), re- 
spectively, where M(T(nS')7r)max is the maximum invari- 
ant mass of the two T{nS)TT combinations. This is used 
to combine T(n5')7r+ and T(nS)'K~ events for visual- 
ization only. Two horizontal bands are evident in the 
T(2S')7r system near 112.6 GeVVc'' and 113.3 GeVVc'', 
where the distortion from straight lines is due to inter- 
ference with other intermediate states, as demonstrated 
below. One-dimensional invariant mass projections for 
events in the T{nS) signal regions are shown in Fig. [2j 
where two peaks are observed in the T{nS)7T system near 
10.61 GeV/c^ and 10.65 GeV/c^. In the following we refer 
to these structures as Z;,(10610) and Zb(10650), respec- 
tively. 

We parameterize the T(55) T(n5)7r+7r~ three- 
body decay amplitude by: 



M = Az, +Az., +Af„ +Af,_ 



where Az-^ and Az2 are amplitudes to account for con- 
tributions from the 2'f,(10610) and Zb(10650), respec- 
tively. Here we assume that the dominant contribu- 
tions come from amplitudes that preserve the orienta- 
tion of the spin of the heavy quarkonium state and, 
thus, both pions in the cascade decay T{5S) — >■ Z{,tt — )■ 
T(nS')7r+7r^ are emitted in an S'-wave with respect to 
the heavy quarkonium system. As demonstrated in 
Ref. 0, angular analyses support this assumption. Con- 
sequently, we parameterize the observed Z;,(10610) and 
^(,(10650) peaks with an S'-wave Breit-Wigner function 

BWis,M,T) = jj^. 

To ac- 



M^-s-iMr ' '^tisre we do not consider 
possible s-dependence of the resonance width. 





■ . ■ (b) 































FIG. 1: Dalitz plots for T(2S')7r+7r" events in the (a) T(2S) 
sidebands; (b) T(2S) signal region. Events to the left of the 
vertical line are excluded. 



with respect to 7r+ and tt transposition. Using isospin 
symmetry, the resulting amplitude is written as 

Az, = az,e'^'^ M^, F^) + BW{s2, Mk,Tk)), 

(2) 

where si = M^[T{jiS)tt+], S2 = AP[T{nS)TT']. The rel- 
ative amplitudes az^-, phases Sz^j rnasses Mk and widths 
Ffc (k = 1, 2) are free parameters. We also include the Afg 
and A amplitudes to account for possible contributions 
in the tt+tt" channel from the /o(980) scalar and /2(1270) 
tensor states, respectively. The inclusion of these two 
states is needed to describe the shape of the M{tt^tt^) 
spectrum around and above M(7r+7r~) = l.OGeV/c^ 
for the T(15')7r+7r~ final state (see Fig. [3]). We use 
a Breit-Wigner function to parameterize the /2(1270) 
and a coupled-channel Breit-Wigner function [10| for the 
/o(980). The mass and width of the /2(1270) state are 
fixed at their world average values ; the mass and the 
coupling constants of the /o(980) state are fixed at val- 
ues determined from the analysis of B'^ 
A/ [/o (980)] = 950 MeV/c2, g^^ . 
Following suggestions in Ref. 
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0.23, OKK 
13, the non-resonant 



amplitude Am is parameterized as A„ 



count for the possibility of T(5S') decay to both Z^tt 
and Z^TT^ , the amplitudes Az-^ and Az2 are symmetrized 



a" 6 2 where S3 ~ M (tt+tt") (53 is not an inde- 
(1) pendent variable and can be expressed via si and S2 but 
we use it here for clarity), a"'', al"^, Sf^ and 62'^ are free 
parameters of the fit. 

The logarithmic likelihood function C is then con- 
structed as 

£ = -2 ^ log(/,igS(si, S2) + (1 - /sig)B(si, S2)), (3) 

where ^(si, S2) is the density of signal events |Af(si,S2)P 
convolved with the detector resolution function, B(si, S2) 
describes the combinatorial background that is consid- 
ered to be constant and /sig is the fraction of signal events 
in the data sample. Both S{si,S2) and B{si,S2) are ef- 
ficiency corrected. 

In the fit to the T(lS')7r+7r^ and T(2S')7r+7r^ samples, 
the amplitudes and phases of all of the components are al- 
lowed to float. However, in the T(3S')7r+7r~ samples the 
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FIG. 2: Comparison of fit results (open histogram) with ex- 
perimental data (points with error bars) for events in the 
T(15') (a,b), r{2S) (c,d), and T(35') (e,f) signal regions. The 
hatched histogram shows the background component. 



available phase space is significantly smaller and contri- 
butions from the /o(980) and /2(1270) channels are not 
well constrained. Since the fit to the T(35')7r+7r~ signal 
is insensitive to the presence of these two components, 
we fix their amplitudes at zero. Due to the very limited 
phase space available in the T{5S) — > T(3S')7r+7r~ decay, 
there is a significant overlap between the two processes 
T{5S) Z+TT- and T(5S') Z^tt+. 

Results of the fits to T(5S') T(nS')7r+7r~ signal 
events are shown in Fig. [21 where one-dimensional pro- 
jections of the data and fits are compared. Numerical 
results are summarized in Table U where the relative nor- 
malization is defined as az^jaz^ and the relative phase 
as (5^2 ~ ■ The combined statistical significance of the 
two peaks exceeds 10 ct for all tested models and for all 
T(niS')7r"''7r^ channels. 

The main source of systematic uncertainties in the 
analysis of T{5S) — > T(?iS')7r+7r^ channels is due to un- 
certainties in the parameterization of the decay ampli- 
tude. We fit the data with modifications of the nom- 
inal model (described in Eq. [T]). In particular, we vary 
the M(7r+7r~) dependence of the non-resonant amplitude 
Anr, include a D-wave component into ^nr, include the 
/o(600) state, etc. The variations in the extracted Zf, 
parameters determined from fits with modified models 
are taken as estimates of the model uncertainties. Other 





10.6 10.7 
M^i,^(n), GeV/c^ 



10.6 10.7 
M^,Ji), GeV/c^ 



FIG. 3: The (a) ht(lP) and (b) ht{2P) yields as a function 
of A/miss (tt) (points with error bars) and results of the fit 
(histogram). 



major sources of systematic error include variation of the 
reconstruction efficiency over the Dalitz plot and uncer- 
tainty in the cm. energy. Systematic effects associated 
with uncertainties in the description of the combinato- 
rial background are found to be negligible. The overall 
systematic errors are quoted in Table HI 

To study the resonant substructure of the T(5S') — >■ 
/if,(mP)7r+7r~ (m = 1,2) decays we measure their yield 
as a function of the ht{lP)n^ invariant mass. The decays 
are reconstructed inclusively using the missing mass of 
the TT+TT" pair, Mmiss(7r+7r~). We fit the Miniss(7r+7r~) 
spectra in bins of /if,(lP)7r* invariant mass, defined as the 
missing mass of the opposite sign pion, Mniiss(7r^). We 
combine the Mniiss(7r"'"7r~) spectra for the corresponding 
-M^miss(7r+) and Mmiss(7r^) bins and we use half of the 
available Mi„iss(7r) range to avoid double counting. 

Selection requirements and the Mijjiss(7r+7r~) fit pro- 
cedure are described in detail in Ref. [3[ . We consider all 
well reconstructed and positively identified tt+tt" pairs 
in the event. Continuum e~^e~ qq {q — w, d, s) 
background is suppressed by a requirement on the ra- 
tio of the second to zeroth Fox- Wolfram moments i?2 < 
0.3 [l3| . The fit function is a sum of peaking compo- 
nents due to dipion transitions and combinatorial back- 
ground. The positions of all peaking components are 
fixed to the values measured in Ref. |3|. In the case of 
the hb{lP) the peaking components include signals from 
T(5S') ^ hbilP) and T(5S') ^ T(2S') transitions, and 
a reflection from the T{3S) — ?> T(IS') transition, where 
the T(3S') is produced inclusively or via initial state ra- 
diation. Since the T{3S) — !• T(IS') reflection is not well 
constrained by the flts, we determine its normalization 
relative to the T(5S') — > T{2S) signal from the exclu- 
sive ^"'"/x~7r+7r~ data for every Miniss(7r) bin. In case of 
the hb{2P) we use a smaller Mniiss(7r+7r~) range than 
in Ref. Mmiss(7r+7r-) < 10.34 GeV/c^, to exclude 
the region of the K'g — >■ tt+tt" reflection. The peak- 
ing components include the T(5iS') — >■ hi,(2P) signal and 
a T(25) — > T(15) reflection. To constrain the normal- 
ization of the T{2S) T(IS') reflection we use exclu- 
7r~ data normalized to the total yield of the 



sive fJ.~^ ^J. 
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TABLE I: Comparison of results on ^(,(10610) and Z(,(10650) parameters obtained from T{5S) 
and T{5S) — >■ hb{mP)-K^-K~ [m = 1,2) analyses. 



T(n5)7r+7r- (n = 1,2,3) 



Final state 



T(15)7r+7r- 



T(25')7r^ 



M[Zt{106W)], MeV/c^ 
r[Z6(10610)], MeV 
M[Zi,(10650)], MeV/c^ 
r[^6(10650)], MeV 
Rel. normalization 
Rel. phase, degrees 



T(3S')7r+7r- hi{lP)n+TT- hb{2P)n 



10605 ± 2"^ 

11 .+4.5+2.1 
^^■^ -3.9-1.2 

10654 ±3tl 



10599+;^+^ 

1 q +10 + 9 

ro _g_7 
106511^11 
19±7l" 

1 f;+0.6+0.4 
^•"-0.4-0.6 

101 +65+74 
-"-''-'--105-109 



10611 ±4±3 
22.3 ± 7.7tl-°o 
10657 ± 6 ± 3 
16.3±9.8l 



+6.0 
2.0 

0.57 ± o.2ii;;j^ 

4 

:9 



58 ± 43: 



10609 ± 2 ± 3 
24.2±3.lt3:g 
10651 ± 2 ± 3 
13.3±3.3l3:o 
0.86±0.1lio'i'o 

+ 17 



10608 ± 2 ± 3 
17.6 ±3.0 ±3.0 
10652 ± 1 ± 2 
8.4 ±2.0 ±2.0 
0.96 ±0.14 



-13 ± 13: 



-9± 19 



+11 



20.9 



0.08 
0.05 



+5.4+2.1 



.39 ± 0.37^^ °^ 



187 



+44+3 



reflection in the inclusive data. Systematic uncertainty 
in the latter number is included in the error propaga- 
tion. The combinatorial background is parameterized by 
a Chebyshev polynomial. We use orders between 6 and 
10 for the hb{lP) [the order decreases monotonically with 
the Afniiss(7r)] and orders between 6 and 8 for the hi,{2P). 

The results for the yield of T{5S) hi,{niP)TT^ tt" 
(to =1,2) decays as a function of the A/miss (tt) arc shown 
in Fig. 131 The distribution for the hh{lP) exhibits a clear 
two-peak structure without a significant non-resonant 
contribution. The distribution for the hb{2P) is consis- 
tent with the above picture, though the available phase- 
space is smaller and uncertainties are larger. We asso- 
ciate the two peaks with the production of the Zf,(10610) 
and Zf,(10650). To fit the Mmiss{7r) distributions we use 
the expression 



\BWi{s, Ml, Ti) ± ae'^BWiis, A'h, T2) + be'^f^. (4) 



Here = Afmiss(7r); the variables Mk, Tk {k 



1,2), 

a, (/), b and tp are free parameters; is a phase-space 
factor, where p (g) is the momentum of the pion origi- 
nating from the T{5S) (Zb) decay measured in the rest 
frame of the corresponding mother particle. The P-wave 
Breit-Wigner amplitude is expressed as BWi{s, M,r) = 

^Mi-^-iM T' ■ Hsre F is the P-wave Blatt-Weisskopf form 



factor F = 



14| . qo is a daughter momentum 



calculated with pole mass of its mother, R ~ 1.6 GcV^ . 
The function (Eq. S]) is convolved with the detector res- 
olution function (cr = 5.2 MeV/c^), integrated over the 
lOMeV/c^ histogram bin and corrected for the recon- 
struction efficiency. The fit results are shown as solid 
histograms in Fig. [3] and are summarized in Table HI We 
find that the non-resonant contribution is consistent with 
zero [significance is 0.3 a both for the hbilP) and hb{2P)] 
in accord with the expectation that it is suppressed due 
to heavy quark spin-flip. In case of the hb{2P) we im- 
prove the stability of the fit by fixing the non-resonant 
amplitude to zero. The C.L. of the fit is 81% (61%) for 
the hb{lP) [hb{2P)]. The default fit hypothesis is favored 
over the phase-space fit hypothesis at the 18 a [6.7 0-] level 
for the hbilP) [^t(2P)]. 



To estimate the systematic uncertainty we vary the 
order of the Chebyshev polynomial in the fits to the 
Afmiss(7r"'"7r~) spectra; to study the effect of finite 
MmissiT^) binning we shift the binning by half bin size; to 
study the model uncertainty in the fits to the Afiniss(7r) 
distributions we remove [add] the non-resonant contribu- 
tion in the hb{lP) [hb{2P)] case; we increase the width 
of the resolution function by 10% to account for possible 
difference between data and MC simulation. The max- 
imum change of parameters for each source is used as 
an estimate of its associated systematic error. We esti- 
mate an additional 1 MeV/c^ uncertainty in mass mea- 
surements based on the difference between the observed 
T{nS) peak positions and their world averages The 
total systematic uncertainty presented in Table U is the 
sum in quadrature of contributions from all sources. The 
significance of the Zh(10610) and .Zb(10650) including 
systematic uncertainties is 16.0 cr [5.6 cr] for the hb{lP) 
[hbi2P)]. 

In conclusion, we have observed two charged 
bottomonium-likc resonances, the ^{,(10610) and 
Zft(10650), with signals in five different decay 
channels, T{nS)TT^ [n = 1,2,3) and hb{mP)Tr^ 
{m ~ 1,2). The parameters of the resonances are 
given in Table |T1 All channels yield consistent re- 
sults. Weighted averages over all five channels give 
M = 10607.2 ± 2.0MeV/c2, F = 18.4 ± 2.4 MeV 
for the Z6(10610) and M = 10652.2 ± l.SMeV/c^, 
F = 11.5 ± 2.2 MeV for the Zh(10650), where statistical 
and systematic errors are added in quadrature. The 
Zh(10610) production rate is similar to that of the 
Zft(10650) for each of the five decay channels. Their 
relative phase is consistent with zero for the final states 
with the T{nS) and consistent with 180 degrees for 
the final states with hb{mP). Production of the Zb's 
saturates the T(5S') — > hb(jnP)-K^ ■k~ transitions and 
accounts for the high inclusive hb{mS) production 
rate reported in Ref. [1]. Analyses of charged pion 
angular distributions favor the = 1+ spin-parity 
assignment for both the Zb(10610) and ^^(10650). Since 
the T(5S') has negative G-parity, the Zb states have 
positive G-parity due to the emission of the pion. 

The minimal quark content of the Zh(10610) and 
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2'(,(10650) is a four quark combination. The mea- 
sured masses of these new states are a few MeV/c^ 
above the thresholds for the open beauty channels B*B 
(10604.6 MeV/c^) and B*B* (10650.2 MeV/c^). This 
suggests a "molecular" nature of these new states, which 
might explain most of their observed properties [15| . The 
preliminary announcement of these results triggered in- 
tensive discussion of other possible interpretations 16- 
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